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O 

^isi ' mation about both the vulcanization transition and the emergent amorphous solid state. Contact 

is made with various physical ingredients that have featured in experimental studies of amorphous 

j^, colloidal and gel systems and in theoretical studies of the glassy state. 
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Certain density correlators, measurable via various experimental techniques, are studied in the 
context of the vulcanization transition. It is shown that these correlators contain essential infor- 
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I. INTRODUCTION AND BASIC INGREDIENTS 



The vulcanization transition (VT) is an equilibrium phase transition from a liquid state of matter to an amorphous 
solid state. It occurs when a sufficient number of permanent random constraints (e.g. chemical crosslinks) — the 
quenched randomness — are introduced to connect the constituents (e.g. macromolecules) whose locations are the 
thermally fluctuating variables. A rather detailed description of the VT has emerged over the past few years, raiiging 
from a mean- field theory of the emerging amorphous solid state ||l|-g| to the critical properties of the VT itself Q . 

The purpose of the present Paper is to investigate the properties of correlators that solely involve the local monomer 
density, and to examine the extent to which such correlators provide access to informative signatures of the VT and 
^ , the emergent amorphous solid state. Along the way, we shall explore the relationship of these density correlators to 
Q ■ various experimental probes, and also discuss their relationship to the diagnostics of "non-ergodic media" studied by 
O I Pusey, Van Megen and collaborators in their work on amorphous states of colloidal and gel systems 0. We note 
that density correlators closely related to the ones we shall be considering also feature in certain recent approaches 
—^ [ to structural glasses g, and we shall touch upon the relationship between our results and those of these recent 
^ . approaches. 

0^ ' Apart from their connections with related studies by other researchers, we are motivated to explore the properties 

C^^ , of density correlators in the context of the VT for the following reason. From the theoretical perspective, the natural 

^' ' collective coordinate from which to view the VT is not the local density; rather it is the amorphous solid order 

parameter, which becomes nonzero as the amorphous solid state is entered and whose correlator-decay properties 

directly mark the onset of amorphous solidification. However, as we shall discuss further below, from the perspective 

^S ' of experiment, the amorphous solid order parameter is rather more elusive than one would like, the most direct way 

'>»; . to measure it being via incoherent quasi-elastic neutron scattering, whereas probes that couple to the density are 

Cd ' more plentiful. For this reason, we wish to examine density correlators in the vicinity of the VT, and the extent to 

H , which they can provide access to both the structure of the amorphous solid state and the long-ranged amorphous 

I ■ solid order-parameter correlations that develop near the VT. 

^ , The approach that we shall adopt to study these density correlators is based upon a minimal model that takes 

into account the integrity of the macromolecules, their thermal position-fiuctuations, the short-range repulsion of 

the constituent monomers, and the permanent random constraints imposed by crosslinking. This minimal model 

has previously been shown to give an accurate picture of the universal properties of the VT, in the sense that its 

. ^ predictions for the mean-field properties of the amorphous solid state have been verified in the computer simulations 

^^ : of Barsky and Plischke P,M, and those of its critical properties that have been elucidated so far (i.e. the percolative 

Vh ' aspects) are in accordance with the predictions of percolation and related field-theoretic approaches |11[| . In order for 

5t 1 our discussion to be concrete and physical, we shall adopt language specific to randomly crosslinked macromolecular 

systems(RCMSs), although our results apply to a broader class of systems. 

Let us now turn to the issue of the order parameter for the VT. This order parameter is crafted to detect and 
diagnose amorphous solidification; it is the following function of n -I- 1 wavevectors {k°, k^, • ■ • , k"}: 

[mY. ds(expik"-Cj(s))^(expiki-Cj(s))^---(expik"-Cj(s))^J, (1.1) 

where N is the total number of macromolecules, Cj(s) (with j — l,...,N and < s < 1) is the position in d- 
dimensional space of the monomer at fractional arclength s along the j*'' macromolecule, (• • ■)y, denotes a thermal 



average for a particular realization x of the quenched disorder (i.e. the crosslinking) , and [• ■ •] represents a suitable 
averaging over realizations of the quenched disorder. As discussed in detail in Ref. IJ], this order parameter does 
indeed detect and diagnose the amorphous solid state. 

Why is the amorphous solid order parameter measurable in neutron scattering experiments? In quasi-elastic neutron 
scattering the incoherent contribution of the scattering cross-section is proportional to 



N ^1 



V / ds exp (iq • Cj (s, 0)) exp ( - iq • c^ (s, t)) ) , 
7=1 -^0 L. 



(1.2) 



where Cj{s,t) is the position of the monome r at time t, the i — > cxo limit of the correlator being proportional (up to 



disorder averaging) to a special case of Eq. (LI), viz., 

— ^ / ds(expiq-Cj(s))x(exp-iq-Cj(s))x . (1.3) 

On the other hand, in several other experimental techniques, such as those discussed below, it is some form of correlator 
involving the local monomer density 

p{r,t)=Y, / dsS{r~c,{s,t)) (1.4) 

that is probed. One frequently-measured correlator is the auto-correlation function of the local density 
(p(x, 0)/9(y, t))x, or equivalently (p(q, 0)p(— q, t))^, where p(q) is the Fourier transform of p{x), i.e., p(q) = 
J (i'^a;p(x) exp(— iq • x) |13]. For example, in neutron scattering experiments this quantity is proportional to the 
coherent part of the quasi-elastic neutron scattering cross-section (see, e.g., Ref. [H, Sec. HIE), and in dynamical light 
scattering experiments, such as those performed on 'non-ergodic' media by Pusey and van Megen wj, this quantity 
is proportional to the intermediate scattering function (also known as the dynamical structure factor) F{k,t). (The 
average over quenched disorder [• • •] in the present work essentially plays the role of the ensemble average (• • •)e 
of Refs. M.) The present theoretical framework is a static equilibrium framework and, as such, is not suitable for 
computing dynamical correlators. However, by using the cluster property (i.e. the fact that the connected correlators 
vanish for t — > oo) we see that the long-time limit of the density-density auto-correlation function is built from the 
equilibrium entity (p(x))^ (p(y))x O'"' equivalently, its Fourier transform {p{q))^ {p{—q))^, an entity that is calculable 
(up to disorder averaging) within our static equilibrium framework. In fact, our approach to the VT is capable of 
calculating precisely this kind of quantity and, therefore, of providing contact with experiments. 

As our results for density correlators are relatively straightforward, we first report the results, deferring the con- 
struction and operation of the necessary theoretical machinery to subsequent sections. Specifically, we find that: 
(i) The usual (i.e. disorder-averaged) density-density correlator [(p(q)p(— q))^] is insensitive to the VT, depending 
only analytically on the constraint density, both at the level of mean-field theory and beyond (i.e. to one-loop order), 
(ii) The density-density correlator involving two thermal averages, [(p(q))x (p(— q))^], is zero in the liquid phase but 
becomes nonzero, continuously, as the system enters the amorphous solid phase. This behavior is a manifestation of 
the freezing-in of random density fluctuations, which is the hallmark of the amorphous solid state. This correlator 
turns out to be proportional the order parameter (at least for weak coupling between the density and the order 
parameter fluctuations). As the order parameter encodes the fraction of localized particles and the distribution of 
localization lengths, this result indicates that these physical diagnostics are accessible via this density-density corre- 
lator. 

(iii) The four-density correlator involving two thermal averages, [(p(k)p(— k));^ (p(q)p(— q))^^^], which can be realized 
as limt_+oo[(p(k, 0) p(— k, 0) p(q, t) p(— q, t))^], becomes long-ranged as the VT is approached from the liquid side. We 
exhibit this phenomenon at the level of mean-field theory. 

(iv) The two density-chanel signatures of the VT given in (ii) and (iii) also provide a means for identifying certain 
critical exponents at the VT, such as the gel-fraction exponent /?, and the correlation-length exponent v. Therefore, 
these density signatures provide another avenue for accessing experimentally the critical exponents of the VT. 

II. FIELD-THEORETIC FORMULATION: MINIMAL MODEL AND COUPLING TO DENSITY FIELD 

We approach the VT by adopting the spirit of the Landau- Wilson scheme for continuous phase transitions. To 
handle the presence of the random constraints we invoke the replica trick and adopt the Deam and Edwards model |l^] 



for the statistics of the quenched randomness (viz. that the statistics of the random constraints are determined by 
the instantaneous correlations of the unconstrained system). Thus, we are led to the need to work with the n — > 
limit of systems of n + 1 replicas. The additional replica, labeled by a = 0, incorporates the constraint distribution. 
With the effective spatial dimensionality thus being determined to be {n + l)d, symmetry considerations lead to the 
following minimal model g, which takes the form of a cubic field theory involving an order parameter field f7(fc) that 
lives on (n + l)-fold replicated d-dimensional space [Oj: 

/ oc - lim n-Mn[Z"] , (2.1a) 

[Z"]^Jv^nexp{-T^^^), (2.1b) 

T^''''{n)=N Y. [-aT+^\k\^)\n{k)\^^Ng ^ f^(fci) ^^(fc^) ^^(^3) <^fc,+fe,+fc3,o • (2.1c) 

feGHRS fci,fe2,fc36HRS 

Here, r is the VT control parameter, which measures the reduced density of random constraints, and the coefficients 
a, b and g depend on the microscopic details of the system. We use the symbol k to denote the replicated wavevector 
{k°, k^, . . . , k"}, and define the extended scalar product k-c by k°-c''+k^-c^ + - • •+k"-c". The symbol k e HRS denotes 
that the summation over replicated wavevectors is restricted to those containing at least two nonzero component- 
vectors k". (We say that this kind of wavevector lies in the higher-replica-sector, i.e., the HRS.) This condition on 
k reflects the fact that no crystalline order (or any other kind of macroscopic inhomogeneity) is present or fluctuates 
critically in the vicinity of the VT. 



We now extend the effective free energy, Eq. ( 2.1c ), by including the field R that is associated with spatial monomer 
density fluctuations. The field R takes as its argument replicated wavevectors having exactly one nonzero component- 
vector ke". (We denote by {e°'}2=o the collection of unit vectors in replicated space, so that, e.g., a generic rephcated 
vector p can be expressed as X]a=oP"^"-) ^^ term the subset of replicated wavevectors having exactly one nonzero 
component-vector the one-replica-sector (IRS) of wave vectors; we term the correspondin g fi elds i?(ke") IRS fields. 
We extend the effective Landau free energy, Eq. ( ^.lc| ), by incorporating the IRS fields |1J], which represent local 
density fluctuations, and add the significant symmetry-allowed cubic term that couples the order-parameter and 
density fields, thus arriving at 

T4n,R)^T^^^{n)+Tl^^R)~^ ^ i?(fcOi?(fc2)f^(fc3)'5fe,+fe,+£3,6, (2.2a) 

ki,k2ems 
^n'^HR)-^ E {r-' + l\H')\RCk)\' + --- ■ (2-2b) 

fcGlRS 

The term T^ (^Rj is the effective free energy for the density fluctuations; in principle, it also includes non-linear 
couplings between the R fields. This effective free energy term already incorporates the effects of the short-range 
repulsion between macromolecules. The parameter r is the correlation length for density fluctuations. (In the context 
of a dense melt, it is simply determined by the monomer density and the effective excluded-volume interaction 
strength Jl^.) The correlation length r remains large and varies analytically (with the constraint density) across the 
VT, and the R field remains a non-ordered field. This is representative of the fact that the disorder- averaged physical 
monomer-density is homogeneous in both the liquid state and the amorphous solid state. 



In addition to the coupling presented in Eq. (2.2a), there is one further term at cubic order that couples the R and 



ft fields, i.e., the vertex QQR consisting of two HRS fields and one IRS field. It can readily be shown by dimensional 



analysis that both this cubic vertex and that given in Eq.( 2.2a) are irrelevant with respect to the fixed points of 
HRS fl theory near d = 6 and, therefore, the critical properties of the VT (i.e. the fixed point structure, the flow 
equation and the critical exponent s) ar e not affected by the coupling to density fluctuations (at least near d = 6). 



Based on the effective free energy ( 2.2a ), our approach is to explore the correlators of the IRS fields (and hence the 
density correlators), taking into account the effects of the VT in the HRS fields by treating what happens in the HRS 
as "input" to be added to the effective free energy of the IRS theory, and working perturbatively (i.e. effectively we 
assume that h is small). 

The reason that we ignore the the cubic coupling flflR, besides its irrelevance in the renormalization-group sense, is 
that it does not contribute to the density correlators that we are interested in (at least to one- loop order). There are 
two points to make in this regard. First, at the mean-field level, the HRS field can be viewed as an external source for 
the IRS field in the cubic coupling fifii?. Due to translational invariance, (r2(fci)ri(A;2))'^" = 0. (We use (• • •)-^" 



to denote a statistical average weighted by T^^^.) Therefore, on average, the term QQR wiU n ot ge nerate a non-zero 



{R) [{■ ■ ■) denotes an average weighted by the repUcated effective free energy presented in Eq. (|2.2a| )]. Second, at the 
one-loop level (and beyond), this term will renormalize the coefficient r~^ (in a singular way) but, as has already 
been shown in App. B of Ref. [^, at least to the one-loop level, there is (in the replica limit) no contribution to the 
density-density correlator coming from HRS critical fluctuations via this kind of vertex. 

In order to help make the physical content of the results that we shall present clear, we pause to give the relationship 
between the physical density correlators and the R correlators: 



(p(k)p(-k)) 



for a — (3, 



lim(i?(ke°)i?(-ke^)) = <^ ^ (2.3) 

(P(k))^ (pi-k))^ ; for a / (3, 



n— fO 



where c denotes that a correlator is connected. (Such connections can be e stab lished by following the replica technique 
that is used in Appendix A of Ref. Q.) On the right hand side of Eq.( |2.3| ), the correlators differ in the locations 
of the thermal averages; on the left hand side they differ in their replica indices, the former being diagonal and the 
latter being off-diagonal in replica space. 

III. FREEZING-IN OF DENSITY FLUCTUATIONS 

Now that we have constructed an extended model containing not only the critical order parameter (i.e. HRS) fields 
but also the noncritical replicated density (i.e. IRS) fields, we proceed to study the effect of critical HRS phenomena 
on the density fields, treating the latter at the tree level. The basic mechanism at work is that the order parameter 
field, which is capable of ordering spontaneously, couples to the density fluctuations via a cubic vertex that is replica- 
off-diagonal as far as the density fields are concerned. A non-zero value of fi = (Jl)'^" , as occurs in the amorphous 
solid state due to spontaneous symmetry breaking in the HRS, contributes replica-off-diagonal terms to the "mass 
matrix" of the i?-field and, hence, leads to the existence of nonzero replica-off-diagonal density-field correlators. 

In order to see this more clearly, we replace n by its expectation value plus fluctuations, i.e., we write 11 == fi + 6fl 
and, hence, arrive at the effective free energy 

T4R,n)^T^^^n)+T'„^^R,n)-^ ^ Rik)R{k)Sn{k)S^^^^^^k^^^, (S.Ia) 

/isGHRS 

^:^^^{R,n) = ±Yl E {ir-' + ^k^)S'^-'-hni-ke}+ke')a-S"-''))Rike'^)R{-kef) + ---. (3.1b) 

( = 5^/3) 

To arrive at this result, we have taken advantage of the facts that both JF„ and f2 are replica (i.e. permutation) 
symmetric, and that f2(— kie"— k2e'^) is macroscopically translational invariant (i.e. it contains a factor of (5ki+k2.o) [pi- 
We now aim to compute the correlator ^_R(ke") i?(— ke'')). By treating /i as a small quantity and expanding 
perturbatively, a direct calculation yields 

(i?(ke") Ri-ke>^)) = (ii'(ke") ii'(-ke'^))^""' + 0{h^). (3.2) 

To obtain the correlator at the tree (in R) lev el, w e neglect the nonlinear self-couplings of R and then invert the 



coefficient matrix of the quadratic term in Eq. (3.1b). Thus, in the n ^ limit, we arrive at 

N 



hm (i?(ke") R{-kef^)) = ^ + ^ NMli~ke^ + kP) „ , . (3-3) 



-*o 



Nhn{~ke^ + ke^ 
2 (r-2 + icP)((r-2 + icfc2) + hTli^ke^ + ke^)) ' 



for a = [i\ 
for a ^ p. 



In the liquid state we have fi = 0, and therefore (_R(ke") i?(— ke^)) = (for a ^ /3). However, in the amorphous 
solid state Ti{kxe^ + k2e2) = q5\,^+i,^,ow{kl + kl) ^ 0, and therefore {R{ke°') R{-ke^)) ^ (for a ^ /3). Here, the 
number q is the gel fraction and the function u;(|/c|2), which decays rapidly with increasing wavevector magnitude on 
the characteristic wavevector scale t^ [pj, encodes the distribution of localization lengths. The simplest setting for 



the density-density correlator eme rges near the VT, where q is small and i(;(|A;p) is negligible unless |fc| < t'^. In this 
regime, by making use of Eq. (|2.3|) we find that 



-| r 0, liquid state; 

(P(k))x (p(-k))^ J = I (^/j^4/2)n(-kei + ke^) = (7V/irV2)r^ wik^r-^"), amorphous solid state. ^^''^^ 

On the other hand, the diagonal correlator [(p(k)p(— k)) ] does not vary with t (and hence varies smoothly with 
the physical constraint density). 

Deeper into the amorphous solid state, the order parameter fl does not decay so rapidly with k and hence the 
quantity [ (p(k)) (p(— k)) 1 should remain appreciable (and thus experimentally accessible) over a wider range of 



k. Now, we expect Eq. (3.3) to remain valid, provided the coupling h is small and both the wavevector dependence 
of h and the finer wavevector dependence of IRS bare correlator are incorporated. (We have omitted the wavevector 
dependence of h so as to simplify our presentation.) Under these circumstances, the wavevector dependence of the 
replica-off-diagonal correlator has the possibility of exhibiting additional features, representative of the ordinary 
density-density correlator [(p(k)p(— k)) ] , superposed on the decaying trend due to the factor w (i.e. due to random 
monomer localization) ||l9(| . 

IV. INHERITED CRITICALITY OF DENSITY CORRELATORS 

In the HRS field theory, the VT is signaled in two ways: (i) via the emergence of a nonzero order parameter, 
and (ii) via the divergence of the correlation length of order-parameter fluctuations. We have already studied the 
replica-off-diagonal density correlator, which is closely related to the order parameter and becomes nonzero upon 
entering the amorphous solid state. We now examine the four-field density correlator [(p(k) /3(— k))^ {p{p) p{—p))^], 
which has the property that it becomes long-ranged at the VT. 




ke^ - pe 



."P 




FIG. 1. Divergent 4-field tree level density correlator. Solid line indicates the bare HRS correlator; wavy lines indicate bare 
IRS correlators. 

We calculate the correlator mentioned in the previous sentence at the tree level (with respect to h vertices) in the 
liquid state and at the VT itself: the Fcynman diagram shown on Fig. n is the only contribution, and gives 

[{p{k)p{-k))^{p{p)p{-p))^] = \im{R{knR{^pef')R{-ke-)R{pe^)) 

oc (/i/7V)2|(i?(ke")i?(-ke"))^""'}'{(i?(pe'3)i?(-pe'^))^"'''}V(ke" -pe^)^(-ke" -t-pe''))^^^^^^ (4.1) 

As anticipated, this density correlator becomes long-ranged at the VT, due to the factor of the HRS order-parameter 
correlator (riri)-^" which, itself, becomes long-ranged at the VT. 

V. DISCUSSION AND CONCLUSIONS 

We have studied density-sector correlators that furnish analogs of the two principal order-parameter signatures of 
the VT: the off-diagonal density correlator [(p(q))x (p(~q))x]' which becomes nonzero as the amorphous solid state 
is entered; and the four-field density correlator [{p(k) p(— k))^ {p{p) p{^P))x]' which becomes long-ranged at the VT. 
We have shown that these density correlators provide useful information about both the emergent amorphous solid 
state and the critical properties of the transition itself. They provide schemes for accessing experimentally the kinds of 
quantities that have been found useful in theoretical investigations of the liquid, critical and amorphous solid states. 



e.g., of vulcanized matter. We are not aware of any explorations of such signatures in the density correlators via 
percolation/gelation approaches to the VT. 

Although, as we have seen, the off-diagonal density correlator [(p(q))x (p(^q))x] ^^ closely related to the order 
parameter of the VT, its non-zero value being induced by a non-zero order parameter, the off-diagonal density 
correlator can be used to diagnose a more general class of systems. By looking at the microscopic definition of 



the order parameter for the VT, Eq. (IT) (especially for the special case in which only two thermal averages are 
involved), we see that it involves only a single summation over monomers and, therefore, the order parameter cannot 
be expressed in terms of local monomer densities. In essence, the order parameter signals the amorphous solid phase 
via its detection of all individual monomers that are localized (and, owing to the explicit crosslinks, distinguishable). 
On the other hand, the off-diagonal density correlator [(p(q))x (p(~q))x] signals the amorphous solid state via the 
detection of the frozen structure of the local density fluctuations. In RCMSs, the localization of monomers invariably 
induces the freezing of this structure, so these two quantities are equivalent. However in systems such as glasses, the 
(local density) structure is, presumably, frozen but each individual particle is able to wander throughout the system, 
given enough time, and therefore all particles retain their indistinguishability Pq |. A quantity essentially identical 
to the off-diagonal correlator considered here has been employed by Mezard and Parisi jg] in the context of their 
theoretical approach to glassy systems. A similar quantity is employed in the diagnosis of the freezing-in of structure 
in colloidal glassy systems and gel systems in the dynamical light scattering experiments of Pusey and Van Megen 
and collaborators [[7|pO[ 

Vulcanized matter certainly differs from glassy systems, inasmuch as it possesses explicit quenched disorder in 
the form of permanent random crosslinks. Because of this, it affords a framework in which the off-diagonal density 
correlator and other density correlators are directly calculable in a well controlled way, both near the transition and 
in the amorphous solid state. Despite the difference with glassy systems, we hope that results such as those presented 
here in the setting of vulcanized matter will not only be useful but also shed some light on the more difficult problem 
of glassy systems. 
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